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Abstract 



Describing configuration spaces of 1-degree-of-freedom (1-dof) linkages in 2D is a difficult 
problem with a long history. We use the Cayley configuration space, a set of intervals of possi- 
ble distance- values for an independent non-edge. We study the following complexity measures 
(associated with the linkage's underlying graph): (a) Cayley size, i.e., the number of intervals, 
(b) Cayley computational complexity of computing the interval endpoints, as a function of the 
number of intervals, (c) Cayley (algebraic) complexity of describing the interval endpoints. 

In both parts of this paper, we restrict ourselves to 1-dof linkages obtained by dropping a bar 
from minimally rigid, tree-decomposable linkages, widely used in engineering and CAD, because 
they are quadratically-radically solvable {QRS, also called ruler-and-compass realizable). 

In Part I of this paper, we give an algorithm to determine the interval endpoints of a Cayley 
configuration space of a 1-dof tree-decomposable linkage, by characterizing the Cartesian real- 
izations corresponding to these endpoints. We then focus on graphs with low Cayley (algebraic) 
complexity. For corresponding generic linkages (as defined here), we show how to find a path 
of continuous motion (provided one exists) between two given realizations, in time linear in a 
natural measure of the length of the path. We show that the number of such paths is at most 
two. In addition, we consider Measures (a) and (b) above for graphs with low Cayley complex- 
ity, and give a natural, minimal realization type, i.e. a minimal set of local orientations, whose 
specification guarantees Cayley size of 1 and 0(|V^P) Cayley computational complexity. As a 
corollary, given two realizations of a linkage with the same minimal realization type, a continuous 
motion path between them is guaranteed, which maintains the same minimal realization type. 
Specifying fewer local orientations results in a superpolynomial blow-up of both the Cayley size 
and computational complexity, provided P is different from NP. 
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1. Introduction 



A linkage {G, I), is a graph G = {V, E) with fixed length bars as edges, i.e. I : E ^M.. A. 
2D realization or configuration G{p) of G is an assignment of points p £ to the vertices 
of G, satisfying the bar lengths in [, i.e., for all edges {u,v) S G, \\pu — Pv\\ = ^(w,!') 
where Pu,Pv S M?. Note that a linkage may or may not have a 2D realization. 

In the following we will use standard and well-known terminology from geometric 
constraint solving and combinatorial rigidity. For detailed background in these fields we 
refer the reader to, for example. Graver et al. (1993) and Sitharam (2005). The degrees 
of freedom (dofs) of a linkage is the minimum number of new bars that need to be added 
in order to ensure that a generic realization of the linkage (with the new bars) is rigid, 
i.e., its bar lengths permit no motions other than the Euclidean or rigid body motions 
(translations or rotations), otherwise the linkage is fiexihle. 

Describing configuration spaces of 1-dof linkages is a difficult problem with a long 
history. In fact, even for rigid linkages, the number of configurations can be exponential 
in the number of vertices and not easy to estimate (Borcca and Streinu, 2004) . Attempts 
to classify configurations according to realization types are given in Sitharam et al. (2006); 
Fudos and Hoffmann (1997), where a realization type uniquely determines a configuration 
of a rigid linkage. For flexible linkages, a well-known early result is Kcmpc (1875) showing 
that an arbitrary algebraic curve can be traced by the motion of a linkage joint. One 
outstanding example is the Peaucellier-Lipkin linkage, which transforms planar rotary 
motion into straight-line motion (Kempe, 1877). For polygonal linkages, recent results on 
the variants of Carpenter's rule problem and pseudotriangulations yield spaces of non- 
crossing configurations and expansive motions (Streinu, 2005; Connelly et al., 2000; Rote 
et al., 2003). Versions of the problem play an important role in Computer- Aidcd-Dcsign 
(CAD), robotics and molecular geometry (Sacks and Joskowicz, 2010; Ying and Iyengar, 
1995; Sitharam, 2005), but only a few results are known beyond individual or specific 
families of linkages (Joan-Arinyo et al., 2003; van der Meiden and Bronsvoort, 2006; Gao 
and Sitharam, 2009; Zhang and Gao, 2006). 

A reasonable way to describe the space of 2D realizations of a 1-dof linkage (G, I) is 
to take a pair of vertices whose distance is not fixed by the bars, i.e., an independent 
non-edge f with G Li f being minimally rigid, and ask for all the possible lengths If 
that the non-edge / can attain (Gao, 2008; Gao and Sitharam, 2005; Sitharam and Gao, 
2008) (i) over all the realizations of (G, /); (ii) over all realizations of (G, I) of a particular 
T -realization type. A T -realization type cr of a 2D realization G{p), where T consists 
of triples of points in p, is a set of local orientations (chirality) at G {+1, —1,0}, each 
denoting the local orientation of a specific triple of points t G T (see Definition 2.3). 

For (i) (resp. (ii)), we call each realizable length If as a (resp. T-oriented) Cayley 
configuration, and the set of all such configurations as the (resp. T-oriented) Cayley 
configuration space of the linkage (G, I) on /, parametrized by the distance //. It is a set 
of disjoint closed intervals on the real line. Aside: the Cayley configuration space over 
/ is actually the projection of the Cayley-Menger semi-algebraic set associated with the 
linkage (G, F) on the Cayley non-edge length parameter / (Cayley, 1841). 



Email addresses: sitharamOcise.ufl.edu (Mccra Sitharam), menghanOcise .uf 1 . edu (Menghan 
Wang), heping. gaoOgmail . com (Heping Gao). 

^ Present address: SDE II at Microsoft, Bellevue, Washington, USA 
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Figure 1. A 1-dof linkage 

For example, the linkage (G, /) in Figure 1 is 1-dof, and adding a bar between any of 
the pairs (wi, 17^+2) would make it rigid. For the example in Figure 2, we choose (vijVs) 
as the non-edge / of (G, /) and give the Cayley eonfiguration space of this linkage over 
/, consisting of the three intervals shown. The figure also shows arbitrarily chosen re- 
alizations corresponding to six Cayley configurations. For the example in Figure 3, we 
take forward realization types (a T-realization type defined in Definition 2.3) a and r 
of realizations (CI) and (C2) in Figure 2 respectively, and show the two corresponding 
forward-oriented Cayley configuration spaces over /. Both oriented Cayley configuration 
spaces consist of a single interval. From these oriented Cayley configuration spaces, we 
show how to track the continuous motion between realizations (C) and (D) as If varies. 

1.1. Questions 

In this paper, we consider various natural questions about the complexity in terms of 
Cayley configuration spaces of a well-known class of 1-dof linkages. 

• (In Part I) How can we obtain an (oriented) Cayley configuration space? How can 
we use an (oriented) Cayley configuration space to obtain canonical continuous paths 
between realizations? Is there a robust measure of complexity of the (oriented) Cayley 
configuration space of a linkage (G, I) over / that depends only on the graph G and 
the non-edge /, and not on the bar lengths 17 

• (In Part II) Is there a robust measure of complexity for (oriented) Cayley configiu-ation 
spaces, that does not even depend on the choice of non-edge /? If so, is there a natural 
characterization and/or efficient algorithmic characterization of graphs G that have 
(oriented) Cayley configuration spaces of low complexity? 

In order to state our contributions more precisely, in the next section, we define mean- 
ingful complexity measures for (oriented) Cayley configuration spaces. 

1.2. Complexity measures for Cayley configuration spaces 

Consider the (oriented) Cayley configuration spaces of 1-dof linkages with underlying 
graph G over a non-edge /. We take the following as measures of complexity: 

(a) Cayley size of G, i.e. the maximum number of intervals in the complete description 
of the (oriented) Cayley configuration space, over all possible linkages (G, I) of G. 

(b) Cayley Com,putational Complexity of G, the maximum overall time complexity of ob- 
taining the complete (oriented) Cayley configuration space, over all possible linkages 
(G, I) of G. It can be regarded as a function of Cayley size. 
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Figure 2. The circled part is the Cayley configuration space of the 1-dof linkage in Figure 1 over 
/ = (vijWa) (demonstrated using our new CayMos software under development (Sitharam and 
Wang, 2012)). Arbitrarily chosen realizations for Cayley configurations: (A) Z/ = 2, (B) If = 3, 
(C) If = 5.5, (D) If = 5.3, (E) If = 8, (F) If = 9. 

(c) Cayley complexity of G, i.e. the maximum algebraic descriptive complexity of each 
endpoint in the (oriented) Cayley configuration space, over all possible linkages (G, I) 
of G. Specifically, if the bar lengths I are in Q, it is desirable if the endpoints are so- 
lutions to a triangularized quadratic system with coefficients in Q (i.e. the endpoints 
belong to an extension field over Q obtained by nested square- roots). Such values 
are called quadratically-radically solvable, or QRS. 
Before investigating these complexity measures, we discuss two desirable requirements, 
completeness and low realization complexity, on the (oriented) Cayley configuration space. 

Completeness means that for each (oriented) Cayley configuration If, there exist only 
finitely many (could be exponential in jGI for the non-oriented case) realizations of (GU 
/, I). This is guaranteed if the linkage (G U /, I) is rigid. 

Low realization complexity means the existence of a linear time algorithm to convert 
from a Cayley configuration to a corresponding Cartesian configuration, with a speci- 
fied realization type (rcsp. from an oriented Cayley configuration If to the corresponding 
Cartesian configuration). 

As an example, the linkage in Figure 1 satisfies both requirements when we choose 
any non-edge / = {vi,Vi+2) as the Cayley parameter, since there exists a simple ruler 
and compass realization of the linkage (G U /, I) from any such /. For linkages with such 
a realization, the coordinate values of realizations are QRS, and we call such linkages 
QRS linkages, and the underlying graphs QRS graphs. 

With these two requirements in mind, we focus on a natural class of 1-dof linkages 
called 1-dof tree-decomposable linkages. The underlying graphs are obtained by drop- 
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Figure 3. Two oriented Caylcy configuration spaces of the 1-dof linkage in Figure 1 (tiiis func- 
tionality is under development in our new CayMos software (Sitharam and Wang, 2012)). Unique 
realizations for oriented Cayley configurations (A1)(A2) If = 5.3, (B1)(B2) If = 5.5 (realizations 
(C1)(C2) from Figure 2), (C) = 5.9, (D) // = 6.3. 

ping an edge from so-called tree-decomposable graphs (formally defined in Section 2). 
Tree-decomposable graphs are minimally rigid and well-studied, for example, in geomet- 
ric constraint solving for CAD. because they are QRS (Fudos and Hoffmann, 1997). 
Conversely, QRS has been shown in Owen and Power (2007) to generically imply tree- 
decomposability in the case of planar graphs, and the implication is strongly conjectured 
for all graphs. Our initial example in Figure 1 is a 1-dof tree-decomposable linkage. 

1.3. Contributions 

Contributions of Part I: In Part I we only consider generic 1-dof, tree-decomposable 
linkages. By a generic linkage, we mean that no bar length is zero, all bars have distinct 
lengths and at most one pair of adjacent bars can be coUinear in any realization. We 
answer the following questions posed in Section 1.1. 

(1) How to obtain an (oriented) Cayley configuration space? 

We answer this question by giving two algorithms for obtaining an (oriented) 
Cayley configuration space. One algorithm works for any 1-dof tree-decomposable 
linkage. The other only works for linkages whose underlying graphs have low Cayley 
complexity - i.e. all interval endpoints in the corresponding (oriented) Cayley con- 
figuration spaces being QRS - using properties of such graphs proved in Part II of 
this paper. 

(2) Given a linkage whose underlying graph has low Cayley complexity, what conditions 
are necessary and sufficient for small Cayley size and, correspondingly, low Cayley 
computational complexity? 

We answer this question by giving a natiiral minimal realization type, i.e. a set 
of local orientations, whose specification guarantees Cayley size of 1 (i.e., the set of 
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realizable distances for the chosen non-edge is a single interval) and Cayley compu- 
tational complexity of 0(|T^p). We also show that specifying fewer local orientations 
than those contained in the minimal realization type results in a superpolynomial 
blow-up of both the Cayley size and computational complexity, provided P is differ- 
ent from NP. 

(3) From (1) and (2) we can immediately answer the following questions for linkages 
whose underlying graphs have low Cayley complexity: given two realizations or two 
Cayley configurations of a linkage, can we determine if there exists a path of contin- 
uous motion between them? How do we obtain such a path if it exists? 

It has been shown in the paper Hidalgo and Joan-Arinyo (2011) that a continuous 
motion path can be obtained by exhaustively searching in a routing graph con- 
structed from what amounts to oriented Cayley configuration spaces (although they 
do not use the Cayley terminology). In our paper, we further show that for generic 
1-dof tree-decomposable linkages with low Cayley complexity, the path between two 
realizations is at most two, and can be directly obtained from the oriented Cayley 
configuration spaces. Provided a path exists, it can be found in time linear in the 
number of interval endpoints of oriented Cayley configuration spaces that the path 
contains. Moreover, when the two realizations have the same minimal realization 
type as in (2) above, it is guaranteed that there exists a path between them staying 
within the same minimal realization type as these two realizations; and the time 
complexity of finding that path is 0(1). 

Contributions of Part II: Part II of the paper is self-contained and answers the follow- 
ing questions from Section 1.1: consider the Cayley configuration space of a 1-dof linkage 
with a underlying graph G over any non-edge /, such that G U / is tree-decomposable. 
Does the Cayley complexity depend on the choice of /? We answer this question in the 
negative. Specifically, we show that if the interval endpoints are QRS for a Cayley con- 
figuration space over some choice of /, then they are QRS for any choice of /. This shows 
robustness of the Cayley complexity measure for 1-dof tree-decomposable graphs. 

Finally, we answer the question: how can we characterize 1-dof tree-decomposable 
graphs G with low Cayley complexity on /, without checking whether every interval 
endpoint in the Cayley configuration space is QRS? We show a surprising result that 
(graph) planarity is equivalent to low Cayley complexity for a natural subclass of 1-dof 
tree-decomposable graphs. While this is a finite forbidden minor graph characterization 
of low Cayley complexity, we provide counterexamples showing impossibility of such finite 
forbidden minor characterizations when the above subclass is enlarged. 

Note. Implementation of the algorithms developed in Part I, Part II and further func- 
tionalities is part of our new CayMos software under development, web- accessible at 
http://www.cise.ufl.edu/~menghan/caymos.html. See screen-shots in Figures 2, 13 and 
14. A different manuscript (Sitharam and Wang, 2012) describes CayMos and includes 
Cayley and Cartesian configuration space analysis and motion analysis of common and 
well-known mechanisms, as well as comparison with other software with related func- 
tionalities. 
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1.4- Organization of Part I 

In Section 2, we give basic definitions related to 1-dof tree-decomposable graphs. 

In Section 3, we associate a special so-called extreme graph with each interval endpoint 
of a (oriented) Caylcy configuration space. We also give the precise definition of low 
Cayley complexity. 

In Section 4 we prove the main result: 0(1) Caylcy size and 0(|yp) Caylcy compu- 
tational complexity when a natural minimal realization type i.e. a minimal set of local 
orientations - is fixed, for a graph with low Cayley complexity. We also show that specify- 
ing fewer local orientations than those contained in the minimal realization type results 
in a supcrpolynomial blow-up of both the Cayley size and computational complexity, 
provided P is different from NP. 

In Section 5, wc give an algorithm to find a path of continuous motion between two 
given realizations from the oriented Caylcy configuration spaces, and show that there are 
at most two such paths. We also show that when those two realizations have the same 
minimal realization type, a path staying within the same minimal realization type can 
always be found in constant time. 

2. Definitions and Basic properties of 1-dof tree-decomposable graphs and 
linkages 

Definition 2.1. A graph G is tree-decomposable if: 

• it is a single edge; or 

• it can be divided into three tree-decomposable components, namely tree-decomposable 

subgraphs d, G2 and G3, such that G = G1LIG2DG3, Gir\G2 = ({u3},0), G2nG3 = 
({1^2}, 0) and Gi n G3 = ({t^i}, 0) where vi, V2 and 113 are three distinct vertices (refer 
to Figure 4(a)). 

A graph G is a 1-dof tree-decomposable graph if there exists a non-edge / such that GU/ 
is tree-decomposable. Such an / is called a base non-edge of G and a base edge of G U /. 

In Figure 4, a tree-decomposable graph is decomposed into three tree-decomposable 

components, and Gi is decomposed into Gn, G12 and G13. 

Note, (i) It is well-known that for a tree-decomposable graph, any decomposition of 
the graph into three disconnected subgraphs by splitting at three distinct vertices is a 
tree-decomposition. 

(ii) A 1-dof tree-decomposable graph G can have many base non-edges. That is, G may 
have non-edges /' ^ / such that both G U / and G U /' are tree-decomposable graphs. 
Wc emphasize that this is different from deleting a different edge h from G U /, which 
also gives an entirely different 1-dof tree-decomposable graph from G. 

Definition 2.2. A 1-dof tree-decomposable graph G can be constructed step by step, 

starting from a given base non-edge /. At the A;*'* construction step, two new maximal 
tree-decomposable subgraphs Ti and T2 sharing a step vertex Vk are appended to the 
previously constructed graph G/(fc— 1): Ti and T2 each has exactly one shared vertex, Uk 
and Wk respectively, with G/(fc — 1) (refer to Figure 5). Vertices Uk, Wk are called the base 
pair of verticesat Step k. We denote this construction step by Vk < {uk & Ti,Wk G T2), or 
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(a) 



(b) 



Figure 4. A graph is tree-decomposable if it can be divided into three tree-decomposable com- 
ponents. 

simply Vk < (ufc, Wfc). We call these maximal tree-decomposable subgraphs T, the clusters 

of G. If a vertex v is shared by m distinct clusters, wc say cdeg{v) = m. 

A tree-decomposable graph can be constructed in a similar way from a given base 
edge. 

For example, refer to Figure 5(a). Construction steps from base non-edge {vo,v'q) are: 
vi < {vo e Ti,v'q e T2), V2 <i {vo e T3,v'o G T4), and V3 < {m e T5,v'o e Te). We have 
cdeg{vQ) = 3, cde5(t;3) = 2. 



Figure 5. (a) A 1-dof tree-decomposable graph with three construction steps from base non-edge 
{vo,v'o). (b) The extreme graph of (a) corresponding to the 3'''^ construction step. 

Note, (i) The decomposition of a 1-dof tree-decomposable graph G into clusters is unique 
(Fudos and Hoffmann, 1996). 

(ii) It follows from the definition that for any construction step Vk < (uk,Wk), the base 
pair of vertices Uk and Wk must lie in two different clusters T„ and in G/(fc — 1). We 
call Tu and the base pair of clusters at Step k. 

A linkage (G, I) with a 1-dof tree-decomposable underlying graph G generically has 
one degree of freedom, and a Cayley configuration space with parameter /. Given If, the 
graph construction of GUf from / clearly yields a QRS realization sequence of (GU/, lUlf) 
from /. Here lUlf : E L) f ^ R. Specifically, for each construction step v < {u, w), given 
Pu, Pw and the lengths l{v,u), l{v,w), p„ can be determined by a corresponding simple 
ruler and compass algebraic solution. The realization may not be unique, since for each 
realization step we may have two possible local orientations. 




(a) 



(b) 
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Definition 2.3. A construction step Vk <\ {ukiWk) of a 1-dof tree-decomposable graph 
G can be associated with a local orientation for the corresponding reahzation step of a 
Unkage {G,l). This local orientation crfc takes a value in {+1,-1,0}, which represents 



the sign of the determinant A = 



Pw Pu 



. A forward realization type is a T-realization 



Pv-Pu 



type, where T is the set of triples of points corresponding to vertices {vk,Uk,Wk), where 
Vk o {ukiWk) represent the construction steps of G from /. A forward- oriented Cayley 
configuration space is a T-oriented Cayley configuration space with respect to forward 
realization type. 

Note. In the following discussion, when referring to oriented Cayley configuration spaces 
of 1-dof tree-decomposable graphs, unless otherwise specified, we always mean the forward- 
oriented Cayley configuration spaces. 

For example, refer to Figure 6 (a) and (b). The graph is 1-dof tree-decomposable 
with base non-edge / = {vo,Vq). The realization step V3 < {v\,V2) has two possible local 
orientations: (a) has 0-3 < 0, while (b) has (73 > by flipping U3 to the other side of 



Figure 6. (a) (b) : Two choices exist for local orientation of realization step < («i , ) from base 
non-edge / = {vo,v'q). (c) (d): realizations for extreme linkage ((5/(3), F"^^). (e) (f): realizations 
for extreme linkage (G/(3),Z"™). 

Remark 1. Since the total number of forward realization types is exponential in the 
number of construction steps (many of forward realization types may be empty, i.e. there 
are no realizations of those types), only by knowing the forward realization type can 
we guarantee linear time complexity for realizing a tree-decomposable linkage, i.e., for 
obtaining a realization from a Cayley configuration of a 1-dof tree-decomposable linkage. 





V2 

(f) 
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3. Extreme graphs and interval endpoints of Cayley configuration space 

We use the notion of extreme graphs and extreme linkages to describe the endpoints 
of the Cayley configuration space of a 1-dof tree-decomposable linkage (G, I) over /. 

Definition 3.1. The fc*'' extreme graph for / of a 1-dof tree-decomposable graph G is 
the graph obtained by adding a new edge ek — {u, w) in G f{k — 1), where u and w are 
the base pair of vertices for the fc*'' construction step Vk < {u, w) from /. We denote this 
extreme graph by Gf{k). We call {u,w) the extreme edge of the extreme graph Gf{k), 
and an extreme non-edge of G. 

Note. It is easy to verify using Laman's theorem (Laman, 1970) that any extreme graph 
of a 1-dof tree-decomposable graph is minimally rigid. 

For example, refer to Figure 5. The 7*'' construction step in (a) is < {v5,vq). Con- 
necting 67 = {v5,vq) in G/(6), we get the extreme graph (7/(7) = G/ (6) U ey in (b). 

Definition 3.2. For a 1-dof tree-decomposable linkage {G,l), the fc*'* extreme linkages 
are (G/(fc),P"") and {G f (k) , l'^'^'^) , where min and max represent the two possible ex- 
treme extensions of I for the extreme edge = (u, w), obtained from triangle inequalities: 
r-'"(efc) := \l{u,Vk)-livk,w)\, ^--'^(efc) := l{u,Vk) +livk,w). 

Note. In realizations of (G/(fc),[™") and (G/(fc), the local orientation at = 0. 

These realizations are sometimes called unyielding realizations. 




Figure 7. When p(w5) and p{vfj) are coincident, l{v5,vs) is not a function of 1{vo,v'q). 

In the following discussion, we denote the Cayley configuration space of a 1-dof tree- 
decomposable linkage (G, F) over / by ^/(G, T), and the oriented Cayley configuration 
space with forward realization type a by $/(G, Z, a). 

The next theorem associates a linkage realization with each endpoint of an (oriented) 
Cayley configuration space. Note that if (G, I) is generic, then for a given Cayley con- 
figuration If and forward realization type a, there exists at most one 2D realization: as 
demonstrated in Figure 7, for any vertex v of G, the point py is not unique only if for 
w's realization step v < (u, w), Pu and Pw are coincident and l{v, u) = l(v, w). 

Theorem 3.3. For a generic 1-dof tree-decomposable linkage (G, /) with base non-edge 
f — (vo,Wq), the following hold: 
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(1) The (oriented) Cayley configuration space over f is a set of closed real intervals or 
empty. 

(2) Any interval endpoint in the (oriented) Cayley configuration space corresponds to 
the length of f in a realization of a,n extreme linkage. 

(3) For any vertex v, is a continuous function of If on each closed interval of the 
oriented Cayley configuration space. Consequently, for any non-edge {u,w), l{u,w) 
is a continuous function of If on each closed interval of the oriented Cayley con- 
figuration space. 

Remark 2. While Theorem 3.3 (2) states that every endpoint of the unoriented Cayley 
configuration space corresponds to an extreme linkage, the converse is not true. As an 
example, refer to Figure 2. Realization (D) is an extreme linkage, but it is not an endpoint 
in the unoriented Cayley configuration space. However we will see in Lemma 5.2 that the 
converse is true for oriented Cayley configuration spaces. 

The proof of Theorem 3.3 follows from elementary algebraic geometry and can be 
considered folklore. However, for completeness, a proof is provided in Appendix A. 

Observation 3.4. (i) Theorem 3.3 gives a straightforward algorithm called ELR (ex- 
treme linkage realization) (described in Appendix B) to obtain the (oriented) Cay- 
ley configuration space for a generic 1-dof tree-decomposable linkage {G,l). The 
algorithm could take time exponential in the Cayley size, 
(ii) Deciding whether the Cayley configuration space over / is non-empty is a NP- 
hard problem, and the Cayley computational complexity is superpolynomial unless 
P = NP. 

Proof, (i) The ELR algorithm works by realizing all the extreme linkages for / con- 
sistent with each forward realization type. More precisely, to obtain an oriented Cayley 
configuration ^•/(G, I, a), we realize all the realizable extreme linkages for / that are con- 
sistent with forward realization type a. Note that since the extreme graphs may not be; 
QRS, realizing each extreme linkage can take time exponential in \ V\ (requiring the solu- 
tion of a general multi- variable system of quadratic equations). Additionally, the overall 
time complexity could be exponential in the actual Cayley size, since many candidate 
endpoints generated during this procedure could finally lead to dead ends, because they 
have no realization for the complete linkage and do not belong to the actual Cayley 
configuration space. The detailed version of this algorithm is in Appendix B. 

(ii) The problem of determining the existence of a realization of a tree-decomposable 
linkage is NP-complete by early results (Saxe, 1979). This problem can be reduced to 
our problem of finding the Cayley configuration space over /, whose decision version is 
whether the Cayley configuration space is non-empty. Clearly, a realization exists if and 
only if the Cayley configuration space is not empty. Therefore, the problem of finding the 
Cayley configuration space over / is NP-hard, and the Cayley computational complexity 
is superpolynomial unless P = NP. □ 

By Theorem 3.3, we can obtain a specific measure of Cayley complexity in terms of 

extreme graphs. Given the promise that there exists a realization of (G, I) corresponding 
to a specific extreme linkage consistent with a given forward realization type, we now 
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consider the following question: what is the algebraic complexity oilf, a potential end- 
point of $/(G, 1)7 I.e., what is the Cayley complexity of G on /? The answer depends on 
whether the extreme graph is QRS or not. If it is not QRS, not only would this adversely 
affect the Cayley complexity i.e. the algebraic complexity of the endpoints, the Cayley 
computational complexity could also be exponential in \V\ (solving general quadratic 
systems) . 

Definition 3.5. A 1-dof tree-decomposable G has low Cayley complexity on base non- 
edge / if and only if all extreme graphs of G for / are tree-decomposable. 

We can easily verify that the 1-dof tree-decomposable graphs in Figure 1, Figure 5 
and Figure 6 all have low Cayley complexity on the given base non-edge. 

4. Finding Cayley configuration spaces for linkages with low Cayley com- 
plexity 

Suppose we are given a 1-dof tree-decomposable graph with low Cayley complexity. 
Any interval endpoint of Cayley configuration spaces of its corresponding linkages can be 
computed essentially using a sequence of solutions of one quadratic equation at a time. We 
may ask: are we also guaranteed to have small Cayley size and low Cayley computational 
complexity even for the oriented Cayley configuration spaces? The answer is yes only 
for a fixed minimal realization type (defined below, which is more restrictive than the 
forward realization type). 

Definition 4.1. For a graph G with low Cayley complexity on /, we can construct each 
extreme graph Gf{k) with the extreme edge as base edge. We call this a reverse con- 
struction. Each realization of (G, I) corresponds to a reverse realization type, a sequence 
of local orientations for each reverse realization step. A minimal realization type consists 
of both a forward realization type and a reverse realization type. 

For example, for the linkage in Figure 6, realizations in (a) and (b) have different 
forward realization typos (thus different minimal realization types). Moreover, since the 
underlying graph has low Cayley complexity on base non-edge / = (fo, Wg), the extreme 
graph G/(3) has reverse construction Vq <\ (wi,W2), ■^'o (^1)^2), where (c)(e) and (d)(f) 
correspond to different reverse realization types (thus different minimal realization types) : 
(c)(e) have vq and v'q on the same side of {v\,V2), while (d)(f) have them on opposite 
sides of (ui, U2). 

Observation 3.4 can be extended to show that in order to guarantee small Cayley size 
and computational complexity, we need to specify the minimal realization type. 

Observation 4.2. For a 1-dof tree-decomposable graph G with low Cayley complexity 
on a non-edge /, when the reverse realization type is specified but the forward realization 
type is unspecified, resp. when the forward realization type is specified but the reverse 
realization type is unspecified, the problem of obtaining the complete description of Cay- 
ley configuration space (decomposition of the Cayley configuration space into a union of 
oriented Cayley configuration spaces) of a linkage (G, I) over / can take time exponential 
in \V\. 
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Proof. When the reverse reahzation type is unspecified, the Cayley size can be exponen- 
tial in \V\. Symmetrically, when the forward realization type is unspecified, the number of 
non-empty oriented Cayley configuration spaces can be exponential in \V\. In Appendix 
C, wc use our initial cixaniplc. the linkage in Figure 1 to demonstrate this exponential 
blow-up in Cayley size and computational complexity. Existing examples (Borcea and 
Streinu, 2004) can also be adapted to show the exponential blow-up. Therefore, in both 
cases, obtaining the decomposition of the Cayley configuration space into a union of 
oriented Cayley configuration spaces takes time exponential in □ 



On the other hand, when the minimal realization type is fixed, for graphs with low 
Cayley complexity, we show below that the Cayley size is 0(1) and the computational 
complexity is 0(|yp). 

4-1- Finding Cayley configuration spaces when the minimal realization type is fixed 

Theorem 4.3. For a 1-dof tree- decomposable graph with low Cayley complexity on a 
non-edge f , if the minimal realization type is fixed, then the Cayley size is 0(1) and the 
Cayley computational complexity is 0{\V\'^). 

Idea of the Proof. Wc first prove the theorem for a special subclass of graphs called 
1-path graphs. Informally, these tree-decomposable graphs have a linearly ordered con- 
struction sequence in a well-defined sense (Definition 4.5). In the case of general tree- 
decomposable graphs, this linear order generalizes to a partial order. The proof for 1-path 
graphs will serve as induction basis for the general case. For the proof of the 1-path case, 
we utilize Lemma 5.9 from Part II (Recursive Structure Lemma) concerning the structure 
of 1-path graphs with low Cayley complexity. Based on this lemma, we obtain a quadri- 
lateral interval mapping ( QIM) algorithm that correctly finds the Cayley configuration 
space (Observation 4.7, 4.8), yielding the proof for the 1-path case (Proposition 4.9). 

Next, we prove the multi-path case, by doing induction on the number of paths of the 
graph. This gives a generalization of the QIM algorithm, which however works only when 
the minimal realization type is fixed (Observation 4.10). From this algorithm, we get the 
proof of the main theorem. 

The structure of the proof is schematically shown in Figure 8. 



Lemma 5.9, Part II 




Observation 4.7 
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(Recursive Structure Lemma) 


— >■ 




QIM(l-path) 


— ^ 


1-path case 





















Observation 4.10 
QIM(multi-path) 



Theorem 4.3 



Figure 8. Structure of proof of Theorem 4.3 
We first give the definition of 1-path graphs. 
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Definition 4.4. Given a 1-dof trcc-dccomposablc graph G, a vertex v is in G"s /asf level 
Lt if: (a) cdeg(v) = 2, i.e. there are exactly two clusters Ti and T2 sharing i;; (b) each of Ti 
and T2 has only one shared vertex with the rest of the graph G' = G \ (Ti UT2) = {«}. 

Definition 4.5. A 1-dof trcc-dccomposablc graph G has a 1-path construction from 
base non-edge / = {vo,v'q) if there is only one shared vertex v in the last level Lt, 
other than possibly ?;o and Uq. As long as there exists a base non-edge permitting 1-path 
construction, wc say the 1-dof tree-decomposable graph G has the 1-path property. See 
Figure 9 for an example. 




Figure 9. (a) A 1-path graph with in the last level Lt- (b) A multi-path graph. 

The following lemma proved in Part II gives a key structural property of 1-path, 1-dof, 
tree-decomposable graphs with low Cayley complexity. 

Lemma 5.9. (Part II) (Recnirsivc Structure Lemma) Let G be a non-trivial 1-path, 1- 
dof, tree-decomposable graph with base non-edge f = {vo,Vq) (by non-trivial we mean that 
G has more than two construction steps). The graph G has low Cayley complexity on f 
if and only if the following conditions hold: 

(1) There are only two vertices vi, V2 directly constructed on {vo,v'q). 

(2) If both vo and v'q are in the last level Lt, G \ {vo,v'q} is a 1-path, 1-dof, tree- 
decomposable graph with low Cayley com,plexity on base non-edge (fi,f2) (see Figure 
10(a)). If only vq is in the last level Lt, let the third construction step be V3 < {ui,u'i). 
Then G\{vo} is a 1-path, 1-dof, tree-decomposable graph with low Cayley complexity on 
base non-edge {ui,u'i) or (u'q,v^) (see Figure 10(b)). 

The proof is given in Part II, Section 5. 

Definition 4.6. A four-cycle of clusters consists of four clusters Ti, T2; T3, T4 such that 
each consecutive pair has exactly one shared vertex. I.e., we have TiHT^ = {f 1}, T^f^T^ = 
{^2}, T3 n T4 = {vz], T4 n Ti = {vi}, where v\, v^, vs, V4 are distinct vertices. Edges 
{vi, W3) and {v2, V4) arc called the diagonal edges of the four-cycle. For any other vertices 
Ui, Uj+i belonging to adjacent clusters Ti, T^+i, {ui,Ui+i) is called a chordal edge of the 
four-cycle. 

Observation 4.7. Given a 1-path tree-decomposable graph G with low Cayley com- 
plexity on / = {uo,Uq), for two consecutive base pairs of vertices that we encounter in 
the construction of G from /, there exists a four-cycle of clusters such that: 

(a) The two consecutive base pairs are the two diagonals of the four-cycle, or: 

(b) One base pair is a diagonal of the four-cycle, while the other base pair is a chordal 
edge of the four-cycle. 
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Proof. We prove by induction on the number of construction steps of G. 

In the base case, G has only one construction step, and the observation is vacuously 
true. 

As the induction hypothesis, assume that the observation is true for all graphs with 
less than k construction steps (fc > 2). For a graph with k construction steps, by the 
Recursive Structure Lemma, there are only two possible cases: (1) exactly 2 step vertices 
vi and W2 are directly constructed on / = {uq^u'q) and both Uq and u'q are in the last 
level Lt\ (2) exactly 2 step vertices vi and V2 arc directly constructed on (mo,Mo) and 
only uo is in the last level Lf. Let 113 < (ui,^^), then eo = {uq,u'q) and ei = {u\,u'i) are 
the first two base pairs of vertices. 

For (1), (wi, u'l) can only be (wi, V2) since both uq and u'q are in Lj, see Figure 10 (a). 
Clearly eo and e\ are the two diagonals of the four-cycle T1T2T4T3. Since G \ {vo,v'q} 
is 1-path, 1-dof, tree-decomposable with low Cayley complexity on {ui,u[) and strictly 
fewer construction steps than G, the induction hypothesis holds for all of its consecutive 
base pairs. That covers all consecutive base pairs of G, proving the induction step for G. 

For (2), let vi < (mq G Ti,u'q e T2), V2 < (wo € T^ju'q e T4). In this case, {ui,u[) may 
not be the other diagonal {vi,V2). However, they must form a chordal edge between T2 
and T4, as shown in Figure 10 (b), since uq is in Lf. G \ {vq} is 1-path, 1-dof, tree- 
decomposable with low Cayley complexity on base non-edge (i) {ui,u[) or (ii) (v^^v'q) 
and strictly fewer construction steps than G, so the induction hypothesis holds for all of 
its consecutive base pairs. Notice that for (ii), {vs,v'q) and {ui,u'i) are two diagonals of 
the same four-cycle T1T2T4T3. That covers all consecutive base pairs of G, proving the 
induction step for G. □ 




Figure 10. For Observation 4.7 

4.1.1. The Quadrilateral Interval Mapping (QIM) algorithm 

From Observation 4.7, we obtain a quadrilateral interval mapping (QIM) algorithm 
for finding the Cayley configuration space over /. We will see the advantage of the 
QIM algorithm over the ELR algorithm (described in Section 3): only by using the QIM 
algorithm can wc attain linear Cayley computational complexity for 1-path graphs of low 
Cayley complexity, when the minimal realization type is fixed. 



15 




Figure 11. (a) an example of ellipse C used in QIM; (b)(c)(d) various cases when mapping with 
ellipse C 

Idea of the QIM algorithm. Consider a quadrilateral with four sides si, S2, S3, S4 

and two diagonals ei, 62- Since the linkage lives in M^, the volume of the tetrahedron 
formed by {si, S2, S3, S4, ei, 62} must equal zero. If we know lengths ?(si), ^(s2), 1(33), 
Z(s4), we can get from the volume equation (a so-called Cayley-Menger determinant) an 
implicit ellipse C relating attainable /(ei) and I (62) values. See Figure 11 (a). So from 
the attainable interval [^'(ei), Z''(ei)] of one diagonal ei, we can obtain the attainable 
intervals of l{e2) by mapping [Z'(ei), ^'^(ei)] on the curve C, and vice versa. 

The curve C used in the algorithm has the following useful properties: (1) one specific 
value ofl{ei) can map to up to 2 distinct corresponding values of l{e2)- Figure 11 (b)(c)(d) 
illustrates several cases in determining the interval for 1(62) from [/'(ei), r(ei)]. This 
mapping can take a set of intervals for ^(ei) into double the number of intervals for Z(e2). 
(2) the overall maximum and minimum points of the curve, Pinin(ei), Pmin(e2), Pmax(ei) 
and Pmax(e2), each corresponds to a change in the minimal realization type. For example, 
the upper left segment of C, from Pinin(e2) to Pmax(ei), corresponds to the realization 
type that the vertices of ei lie on different sides of the line specified by 62, and the vertices 
of 62 lie on the same side of the line specified by ei. 

Algorithm 1. QIM: given a linkage (G, I), where G is a 1-path graph with low Cayley 
complexity on base non-edge /, we start from the interval [^™'"(efe), Z™^''(efe)] obtained 
by triangle inequality from extreme linkages of the last extreme graph G/(fc), and map 
back to obtain a set Sq of intervals for If. 

Sk ^ {r'"(efe),r'"^''(efe)]} 
for i = fc — 1 to 

[Map the interval of l{ei+i) back to l{ei)\ 
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if Ci+i and satisfy the requirement in Observation 4.7 (a) 
Obtain Si from Si+i using mapping with the ehipse C 
[ei+i and Cj are two diagonals of a quadrilateral] 
else [si+i and satisfy the requirement in Observation 4.7 (b)] 

if ej+i and Cj connect the same pair of adjacent clusters in the four-cycle 
Obtain Si from Si+i by applying the law of cosines 
[see Remark 3(i)] 

else [cj+i and ej connect different pairs of adjacent clusters in the four-cycle] 
Obtain Si from Si+i by applying the law of cosines and mapping with C 
[see Remark 3(ii)] 
return 5*0 

Remark 3. (i) Refer to (m2, u'2) and [ui^u'^) in Figure 10 (a). Consider the two triangles 
l\v^uiu'i and Av^U2u'2. Since and Tg are fixed clusters, the lengths of triangle edges 
{vs,u\), (t;3,ui), (v3,U2) and {vi,u'2) are fixed. Moreover, if we know one of the two 
angles, /.uiv^u'i and we can easily obtain the other. So from a value of u'l), 

by the law of cosines, we can obtain /.uiv^u'i and thus ZM2W3M2, from which we can get a 
unique corresponding value of l{u2, u'2). Symmetrically, each value of Z('U2, W2) corresponds 
to a unique value of l(ui,u'i). 

(ii) Refer to {u2,U2) and {ux,u'i) in Figure 10 (b). By (i), there is a one-to-one cor- 
respondence between Z(i)3,Uo) and Z(u2,W2) with constant time cost. Notice that (^3,^0) 
and {ui,u'i) are the two diagonals of quadrilateral U3U1U0U1. Therefore we can map to 
l{ui,u{). 

Observation 4.8. For a 1-path tree-decomposable linkage (G, I) where G has low Cayley 
complexity on /, the set Sq of intervals for If generated by the QIM algorithm is exactly 
the Cayley configuration space over /. 

Proof. We prove by induction on the number of construction steps of G. 

In the base case, there is only one construction step, and the Observation is vacuously 
true. 

As the induction hypothesis, we assume that the algorithm correctly generates the 
Cayley configuration space for linkages with less that k construction steps. For a graph 
G with k construction steps, by the Recursive Structure Lemma, there are only two 
construction steps directly based on / = (mq, Wq), and G" = G\{uo, u'q} (or G' = G\{uo}) 
is a 1-path, 1-dof, tree-decomposable graph with low Cayley complexity on {u\, u'{), and 
less than k construction steps. 

By induction hypothesis, (G' ,1) is realizable if and only if l{ui,Ui) is in the set Si of 
intervals generated by applying QIM on {G',l). For QIM on {G,l), we do an additional 
mapping to get Sq from ^i. When If belongs to an interval of Sq, the first four-cycle 
vou[v'qUi is clearly realizable; moreover, l{ui,Ui) is in Si since So is generated by mapping 
from ^i, so (G', I) is also realizable. Thus for all If values in 5o, (G, I) is realizable. On the 
other hand, when If is not in any interval of Sq, either the first four-cycle is not realizable, 
or {G',l) is not realizable, so all realizable If values are contained in 6*0. Therefore So is 
the Cayley configuration space of (G, I) over /. □ 

We give two examples to demonstrate how the QIM algorithm works. 
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Example 1: 

To obtain the Cayley configuration space on / = {vq, v'q) for G in Figure 9(a): 

• Step 1: Obtain the interval of l{v4,V2) in by triangle inequality; 

• Step 2: In quadrilateral v'qV2V3V4, obtain the interval of 1{v'q,V3) from the interval of 

• Step 3: Similarly, in quadrilateral v'qViV3V2, we have 1{vq,V3) — >■ l{vi,V2); 

• Step 4: In quadrilateral Woi'2i'gDi, we have l{vi,V2) 1{vq,Vq). 

Example 2: 

To obtain the Cayley configuration space on / = {vq, v'q) for G in Figure 12: 

• Step 1: Obtain the interval of l{u2,V3) in AW4W3M2 by triangle inequality; 

• Step 2: In four-cycle T^TsT^Tl, we have l{u2,V3) l{ui,u[); 

• Step 3: In four-cycle T^T{TiT2, we have l{ui,u[) 1{vq,v'q). 



The worst-case time complexity of the QIM algorithm is exponential in the number 
of construction steps, since each mapping could possibly double the number of intervals. 
Nevertheless, when the minimal realization type is fixed, using the QIM algorithm, we 
can obtain the Cayley configuration space of a 1-path graph in linear time. 

Proposition 4.9. For a 1-path, 1-dof, tree- decomposable graph G with low Cayley com- 
plexity on f, if the minimal realization type is fixed, for any linkage {G,l), the Cayley 
configuration space over f contains at most one interval, and can he found in 0{\V\) 
time. 

Proof. When the minimal realization type is fixed, when applying QIM to any linkage 
with underlying graph G, we are restricted to a monotonic segment of the ellipse C. 
We start from a single interval for the last extreme edge of G and map back to /. 
Since we are always mapping using monotonic segment of C, the image is at most one 
interval. Therefore, If , as well as each of the distances between base pairs of vertices, 
takes values in at most one attainable interval (the interval depends on the base pair). 
So each construction step takes constant time. Since the number of construction steps is 




Figure 12. QIM for a 1-path graph. See Example 2. 



0{\V\), <^f{G,l) can be computed in 0{\V\) time. □ 
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Next, wc will look at the general ease, where a graph can have more than one path, 
namely, more than one vertex at the last level Lt (other than the endpoints of the given 
base non-edge). We say that each such last level vertex corresponds to a path. 

Observation 4.10. For a 1-dof tree-decomposable graph G with low Cayley complexity 
on / = (i)o,Uo), if the minimal realization type is fixed, we can use the QIM algorithm 
to correctly find the Cayley configuration space over / for any linkage {G,l), and this 
Cayley configuration space is a single interval. 

Proof. Consider a linkage {G, I) with k + 1 paths. For any vertex Vi in the last level Lt 
different from vq and v'q, there exists a sequence of construction steps from f to Vi, such 
that Vi depends on all these construction steps. We call this construction sequence the 
i*^ path Pi of G, which is in fact a 1-path subgraph. By Proposition 4.9, for each pi, we 
can apply the QIM algorithm for 1-path linkages to find ^f{pi, I) over /, which is a single 
interval /j. Then we take the intersection I = f] li, which is also a single interval. 

l<i<k 

By induction hypothesis, each path {pi,l) is realizable if and only if If is in li, so if 
If ^ I, the linkage is not realizable. On the other hand, since the minimal realization type 

is fixed, when ?j G /, all paths are simultaneously realizable (each value of If corresponds 
to a unique realization of the entire linkage). Therefore, / is the Cayley configuration 
space of (G, I) over /. □ 

Finally, we can prove the main theorem. 

Proof, (of Theorem 4.3) By Observation 4.10, for any linkage with underlying graph 

G, the Cayley configuration space over / contains at most a single interval. The QIM 
algorithm takes 0(|V^|) time on each path, and there can be at most O(IV'I) paths, so 
the Cayley computational complexity is 0(|y|^). □ 

Remark 4. This result shows the advantage of the QIM algorithm over the ELR algo- 
rithm we mentioned in Section 3. First, realizing each extreme linkage takes 0(|y|) time, 
therefore realizing all 0(|y |) extreme linkages from / takes 0(|y p) time. Moreover, note 
that when the reverse realization type is fixed, an interval endpoint can also arise from 
a change in reverse realization type. This is not an extreme linkage for the given base 
non-edge, but an extreme linkage for an extreme non-edge. So the ELR algorithm must 
consider extreme linkages for all these 0(|V^|) possible base non-edges, and the overall 
time complexity is 0(|yp). 

An implementation of the algorithms discussed in this paper for finding Cayley con- 
figuration spaces for linkages with low Cayley complexity can be found in the software 
CayMos under development (Sitharam and Wang, 2012). A screen-shot is included in 

Figure 13. 

Why does the QIM algorithm fails as is, if the minimal realization type is not fixed? 
Notice that when the minimal realization type is not fixed, we cannot simply take in- 
tersection of intervals as in proof of 4.10, since realizations requiring different minimal 
realization types can generate the same //. For example, refer to Figure 9 (b). There are 
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trace a vertex 



find motion path 
[Sj tietween two realiiations 
O between Cayie> configurations 




Figure 13. Finding Cayley configuration spaces using our new CayMos software under develop- 
ment (Sitharam and Wang, 2012). 1: labels showing information of the input linkage; 2: options 
for finding the oriented / non-oriented Cayley configuration space, and a spinner for navigating 
the Cayley configuration space (using length of the base non-edge); 3: intervals of the Cayley 
configuration space, where the black point shows the current Cayley configuration; 4: button for 
changing the solution type of the linkage. 

two paths, po and pi, corresponding to vertices and vq in the last level respectively. 
For certain I, pi may require vi and V4 to be on the same side of (v^^Vq) such that vq 
can be realized. At the same time, po may require vi and 1)4 to be on different sides of 
(w3,Uq) such that V5 can be realized. These two different realization types can generate 
the same If. If we just take intersection of intervals of values for If, even the intersection 
is non-empty, the linkage may still be unrealizable since po and pi require conflicting 
realization types. An alternative strategy would be to follow the construction steps, and 
perform 0(|1^|) mappings altogether. 

The above discussion leads to the following conjecture: 

Conjecture 4.11. Even when the minimal solution type is not fixed, the QIM algorithm 
can be adapted to work for 1-dof tree-decomposable linkages whose underlying graphs have 
low Cayley complexity When the minimal solution type is fixed, this adapted algorithm 
can be used to obtain the Cayley configuration space in 0{\V\) time. 



5. Finding a patii of continuous motion between two realizations 

As we will show below, we can easily find a path of continuous motion between two 
given realizations, given the oriented Cayley configuration spaces of a linkage whose 
underlying graph has low Cayley complexity. 
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Theorem 5.1. For a generic linkage [G.l) where G has low Cayley complexity on f, 
(i) Given two realizations, there exist at most two paths of continuous motion between 
them, and the time complexity of finding such a path (provided one exists) is linear 
in the number of interval endpoints of oriented Cayley configuration spaces that the 
path contains. 

(a) Given two realizations with the same minimal realization type, there exists a unique 
path of continuous motion between them, staying within the same minimal realiza- 
tion type, and the time complexity of finding that path is 0(1). 

Idea of the proof. The genericity is very important here. We analyze the continuous 

motion of the hnkagc in Lemma 5.2 and 5.3, proving that the reahzation type can only 
be changed via extreme linkages at interval endpoints of oriented Cayley configuration 
spaces, and (i) directly follows. Then we apply Theorem 4.3 to get (ii). 

As noted in Remark 2, an extreme linkage realization may be an internal point of an 

interval in the unoriented Cayley configuration space. We now show that every extreme 
linkage realization corresponds to an interval endpoint of the oriented Cayley configura- 
tion space. 

Lemma 5.2. (i) Between two points lying in the same interval of an oriented Gayley 
configuration space, there always exists a unique path of continuous motion staying within 

that interval. 

(ii) No extreme linkage realization can correspond to an internal point of an interval in 
any oriented Cayley configuration space. 

Proof, (i) Let If = h and If = I2 he in the same interval la in oriented Cayley con- 
figuration space ^f{G,l,a-), and the corresponding realizations be G{pi) and G{p2) re- 
spectively. By definition of oriented Cayley configuration space, there exists a ruler and 
compass realization with realization type a for every Cayley configuration in 1^ between 
li and I2. One can easily verify that this series of realizations gives a path of continuous 
motion between G{pi) and G(p2). 

(ii) Assume some extreme graph Gj{k) has a corresponding extreme linkage realization 
G(p), with /'s length I j = Iq being an internal point of an interval in oriented Cayley 
configuration space ^ f{G,l,a). Since the length l^ of the extreme edge e in G{p) is the 
max;imum (resp. minimum) possible value for e's length l{e), it follows that for all small 
enough e, the oriented Cayley configurations If = Iq — e and If = Iq + e m both 
correspond to realizations with l{e) < I* (resp. ?(e) > /*). 

By the proof of Theorem 4.3, l{e) would be a monotonic function oilf if the minimal 
realization type is unchanged. Since l{e) increases before If reaches Iq and decreases after 
If reaches Iq, the minimal realization type must change at Zq. Since the forward realization 
type is unchanged, we can conclude that the reverse realization type must change at Iq. 
Namely, there are two pairs of collinear bars in G(p), one corresponds to the k*^ forward 
construction step, the other corresponds to the change in the reverse realization type. 
However this violates our assumption of genericity of linkages that at most two bars can 
be collinear in any realization. □ 
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Lemma 5.3. During continuous motion, a linkage whose underlying graph has low Cay- 
ley complexity can only change forward realization type via endpoints of oriented Cayley 
configuration spaces, and only one entry of the forward realization type is switched by 
such a change. 

Proof. Suppose the linkage changes forward realization type at If = Iq during the con- 
tinuous motion. In other words, the linkage realization has forward realization type a 
immediately before reaching and has forward realization type t immediately after 
reaching Iq, where a ^ t. To guarantee continuous motion, the forward realization type 
at ^0 must be compatible with both a and r. Due to the genericity assumption on link- 
ages, this is possible only if a and r differ at exactly one entry which is at Iq- Therefore 
the realization at lo is a realization of an extreme linkage. By Lemma 5.2 (ii), Iq must be 
an interval endpoint of both I, a) and T, r). □ 

Proof, (of Theorem 5.1) For (i), the following algorithm finds a path of continuous 
motion between two realizations in time linear in the number of endpoints along that 
path: 

From the starting realization G{ps) with forward realization type a, we take the ori- 
ented Cayley configuration ^f{G,l,<T) and find the interval that G{ps) is in. By 
Lemma 5.2 (i), inside I„ there always exists a path of continuous motion. Take one end- 
point If = lo of /(J. In the corresponding realization G{po), exactly one entry of the 
forward realization type is 0. By Lemma 5.3, the next immediately reachable oriented 
Cayley configuration space ^/(G, /,t) is uniquely determined, since that entry should 
have the opposite sign in r. Since intervals in an oriented Cayley configuration space 
are all disjoint, there is at most one interval in i>/(G,r, r) with Iq as an endpoint. So 
we can find at most one interval Ir immediately reachable. We repeat the process from 
the other endpoint of I^- until we reach the interval containing the target realization, or 
we may arrive at a deadend, in that case we start again from the other endpoint of !„■ 
Each endpoint encountered leads to at most one next immediately reachable interval, so 
backtracking is never necessary. Therefore the time complexity is linear in the number 
of endpoints along the path we found. Since both endpoints of I a- could potentially lead 
to the target realization, there are at most two paths between two given realizations. 

For (ii), by Theorem 4.3, the two realizations lie in a single interval la- in the corre- 
sponding oriented Cayley configuration space. By Lemma 5.2 (i), within a single interval 
of an oriented Cayley configuration space, there is always a continuous motion path. 
Therefore, there always exists a uniqiie path maintaining the minimal realization type 
between the two realizations within I^. The time complexity of finding this path is 
0(1). □ 

Figure 3 gives an example running of the algorithm described in the proof above, 
finding a path of continuous motion from realization (Bl) with realization type a to 
(B2) with realization type r. We start from the interval containing (Bl) in $/(G, I, a) 
and take one endpoint of which corresponds to extreme linkage realization (Al). 
Taking the entry of a which is in (Al), and reversing its sign in a, we get the next 
realization type r. Now we go from (Al) to (A2), which is essentially the same realization 
but contained in $/(G, F, t), and the immediately reachable interval Ir with (A2) as an 
endpoint realization is uniquely determined. Since Ir contains the target realization (B2), 
a continuous path is successfully found. 



22 



Remark 5. 1. For Theorem 5.1 (ii), a second path of continuous motion could exist, if 
we remove the requirement that the path should stay within the same minimal realiza- 
tion type. In that path, the linkage leaves la from the endpoint closer to the starting 

realization than the target realization, takes various minimal realization types along the 
path, and reaches the target realization via the other endpoint of 1^. By Theorem 5.1 
(i) , the time complexity of finding this path is linear in the number of interval endpoints 
contained in it. 

2. Not much improvement of the length of the paths is possible beyond Theorem 5.1 
(i), if the two realizations have either the same forward realization type or the same 
reverse realization type. When the two realizations have the same forward realization 
type, the two realizations may belong to different intervals of the corresponding oriented 
Cayley configuration space. When the two realizations have the same reverse realization 
type, the two realizations may belong to different oriented Cayley configuration spaces. 
In both cases, it is not guaranteed that a path of continuous motion exists, and even if 
such a path exists, the number of interval endpoints contained in it is hard to determine. 

Corollary 5.4. To obtain a continuous path between two Cayley configurations where 
their forward realization types are unspecified, we run the algorithm given by Theorem 
5.1 (i) for each candidate forward realization type of the starting Cayley configuration 
and each candidate forward realization type of the target Cayley configuration. 

Note. For each pair of starting and target realizations there are at most two paths, but 
the number of such pairs could be exponential in the size of the linkage. 

An implementation of the algorithms discussed above for finding a continuous motion 
path between two linkage realizations / two Cayley configurations can be accessed in 
the software CayMos under development (Sitharam and Wang, 2012). A screen-shot is 
included in Figure 14. 

6. Conclusion 

In Part I of this paper, we investigated the structure of Cayley configuration spaces for 
1-dof tree-decomposable linkages, introduced the complexity measures associated with 
the underlying graphs, and formally gave algorithms to obtain such Cayley configura- 
tion spaces. For linkages whose underlying graphs have low Cayley complexity, we gave 
sufficient and necessary condition for small Cayley size and low Cayley computational 
complexity. We also gave an efficient algorithm to find a path of continuous motion for 
such linkages. 

Implementation of the algorithms developed in Part I, Part II and further func- 
tionalities is part of our new CayMos software under development, web- accessible at 

http://www.cise.ufl.cdu/"'mcnghan/caymos.html. Sec screen-shots in Figures 2, 13 and 
14. A different manuscript (Sitharam and Wang, 2012) describes CayMos and includes 
Cayley and Cartesian configuration space analysis and motion analysis of common and 
well-known mechanisms, as well as comparison with other software with related func- 
tionalities. 
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Figure 14. Preview of finding a continuous motion path using our new CayMos software under 
development (Sitharam and Wang, 2012). 1: options for finding a continuous motion path be- 
tween two realizations / two Cayley configurations, and buttons for picking the starting and 
ending realization / Cayley configuration; 2: a spinner for traclcing the motion of the linkage 
along the continuous motion path; 3: the linkage shown in the middle tracks the continuous mo- 
tion path, and the starting / ending realizations are shown as thumbnails in the upper corners; 
4: oriented Cayley configuration spaces encountered along the path; 5: button for tracking the 
trajectories of the linkage's vertices; 6: the trajectories of vertex and ug. 
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A. Proof of Theorem 3.3 

Proof. We first prove for a fixed forward realization type a, that the theorem holds 
for the oriented Cayley configuration space, by induction on the number of construction 
steps from /. 

In the base case, G has only one construction step vi < [vq E Ti. E T^)- The 
distances wq) and 1{vi,v'q) are fixed by clusters Ti and T2 respectively, and by triangle 
inequality, $/(G,r,CT) is a single closed interval [\l{vi,V()) — l{vi,v'f^)\A{vi,va)+l{vi,v'Q)\. 
Clearly, (1) and (2) hold. For (3), without loss of generality. let ]3(uo) be the origin, 
p{vq) lie on the x-axis, and p(vi) = {xi, yi) has positive y-coordinate. Let Ri = 1{vi^vq), 
i?2 = 1{'Vi,Vq) and i?3 = ?(vo,fo) = If - In AwoWqWi, we have 

lii -f- Itrt — Itn 

^ + R2 + R3){Ri + R2- R3){Ri -R2 + Ri){-Ri + R2 + R3) 

2i?3 

where cti is the entry corresponding to the first construction step in the forward realiza- 
tion type a. 

Since the linkage is generic, i?3, namely If, cannot be 0, so both xi and yi are contin- 
uous functions of I/. Moreover, since internal realizations of both Ti and T2 are uniquely 
specified, the coordinates of all other vertices in Ti and T2 are continuous functions of 
coordinates of vi, vq and v'q, thus continuous functions of 

As induction hypothesis, assume that the theorem holds for linkages whose underlying 
graph Gj{k — l) has fc — 1 construction steps. Consider a graph Gf{k) with k construction 
steps, obtained by adding one more construction step Vk < (wfe € Ti,Wk G T2) to G f(k — 
1). For any linkage {Gf{k),l), according to Statement (3) of the induction hypothesis, 
l{uk,Wk) is a continuous function of If, say l{uk,Wk) = g{lf)- By triangle inequality, 
l{uk,Wk) is restricted to the interval [min, max] where min — \l{uj^,Vk) — l{wk,vi^)\ and 
max = l{uk, Vk) + l{wk, vi~). This restriction may create new candidate interval endpoints 
in ^f{Gf{k), I, a), namely g~^{mm) and ^"^(max). as shown in Figure A.l. A candidate 
endpoint is actually a new interval endpoint, only if its corresponding extreme linkage 
realization p(G/ (A;), r™, a) (resp. p(G/(/c), P"'^'^, ct)) does exist. So (1) and (2) also hold 
for {Gf{k)J). 

To prove (3), take any vertex v in G f{k). By induction hypothesis, if f ^ (riUT2), p{v) 
is a continuous function of Z/. For v € (Ti UT2), we first consider p{vk). For convenience, 
first rotate and translate the coordinate system so that p{uk) is at the origin, p{wk) 
is on the x-axis, and p{vf.) = {xk,yk) have positive y-coordinate. Let Ri = l{vk,Uk), 
i?2 = l{vkiWk) and R^ = l{uk,Wk)- In AukWkVk, we have 
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h 



h 



Figure A.l. For Theorem 3.3. New constraint on l(uk,Wk) creating interval endpoints in 
'^f{Gf{k)J)- •■■ 5"^(min); o: ^-^(max). 

^ + B2 + R3){Ri + i?2 - R3){Ri - R2 + Rz){-Ri + R2 + R3) 

where is the entry corresponding to the fc*'' construction step in the forward reahzation 
type cr. 

Since the hnkage is generic, R3 > 0. So both Xk and y/c are continuous functions of 
Moreover, since internal reahzations for both Ti and T2 are specified, the coordinates of 
any other v € (T1UT2) are continuous functions oip{uk), p{vk) aiidp{wk), thus continuous 
functions oflf. Consequently, for any non-edge (u, w), l{u, w) is also a continuous function 
oilf- Note that this continuity is not affected even if we transform back to the original 
coordinate system. 

For the complete Cayley configuration space, (1) and (2) still hold since $/(G, is 
just the union of oriented Cayley configuration spaces over all possible forward realization 
types. □ 



B. Finding Cayley configuration spaces by realizing all extreme linkages 
(ELR) 

Given a linkage (G, /) and a forward realization type a, in the ELR algorithm, we 
use a set Ig- to store candidate intervals of the oriented Cayley configuration space, 
which is initially the entire M^. For each realization step i, we update Xo- by considering 
restrictions on If from all extreme linkages realizations of Gf{i) with forward realization 
type (T. After we have done this for every realization step, 1^ is the oriented Cayley 
configuration space. 

Algorithm 2. ELR 

2ct ^ (—00, +00) 

for i = 1 to k do [k is the number of G 's construction steps] 
^ [set of candidate interval endpoints] 
for every extreme linkage realization p of Gf{i) with forward realization type a 
if {GU f,lU If) is realizable 
add I f value of p to 5 
for each candidate endpoint Iq in S 
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UPDATE(Ict, ^o) [see following discussion] 
return 



When updating To-, we need to notice that not every candidate Cayley configuration 
^0 in S actually creates new restriction on I<j- Recall from the proof of Theorem 3.3 
that a realization step vt < {ui,Wi) restricts l{ui,Wi) = g{lf) in [min, max]. As shown in 
Figure B.l, there are three possible cases for a candidate configuration Iq: (a) both the 
left and the right neighborhood of Iq fall into 1^; (b) the left neighborhood of Iq falls 
into but the right does not, and symmetrically, the right neighborhood falls into 
but the left does not; (c) neither the left nor the right neighborhood of Iq falls into To-, 
meaning that If itself is the only realization in the neighborhood. In (b), creates an 
new endpoint in I^- In (c), lo creates an isolated point in In (a), Zq does not create 
any interval endpoint in ■ 

^l{Uk,Wk) nl{uk,Wk) tfl(Uk,Wk) 

g{lf) 



l{Uk,Wk) 




Figure B.l. min = \l{vk,Uk) — l{vk,Wk)\, max = l{vk,Uk) + l{vk,Wk)- (a) candidate If values o 
that are internal points of some interval in not endpoints; (b) candidate If values • that are 
new interval endpoints in creating intervals 7i and I2; (c) candidate If value o that creates 
an isolated point in X^. 

So in the UPDATE procedure, for each candidate Cayley configuration Iq, we check if 
there is any realization, with If value between Iq and the immediately preceding (resp. 
immediately succeeding) candidate interval endpoint in 



/min /max 



Algorithm 3. UPDATE(X^ 

prev ma,x{l Gla\l < Iq}, next <r- min{Z £1^1' > lo} 
p <r- {prev + Zo)/2, n <— {next + Iq)/2 

P true if p has corresponding realization, false otherwise 
N true if n has corresponding realization, false otherwise 
if exactly one of P and TV is true 

add Iq as an endpoint in 
elseif both P and N are false 

add lo as an (isolated) endpoint to 



To obtain the complete Cayley configuration space, we just take ^f{G, I) = (JX^. 



C. Exponential blow-up of Cayley size and Cayley complexity when the 
minimal realization type is not specified 

We provide in this section an example, that the Cayley size of a graph with low Cayley 
complexity can be exponential in \V\, even when a forward realization type is specified. 
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One can use a symmetric example to show that if we just fix the reverse reahzation type, 
there can be exponentially many non-empty oriented Cayley configuration spaces. 

Observation C.l. The Cayley size of a graph with low Cayley complexity can be 
exponential in the number of construction steps if we only fix the forward realization 
type. 

Proof. We give an example of a 1-dof tree-decomposable linkage {G,l) with a fixed 
forward realization type a and low Cayley complexity, which has exponential Cayley 
size. See Figure C.l. The base non-edge is / = /i = (i^i, fa). For convenience we slightly 
abuse our notation to let the construction step number start from 0, so the 0*'' and 1*** 
construction steps are V4 < {vi,V3) and V2 < (t^i, fa) respectively. They form the outermost 
quadrilateral Qi = V4V3V2V1. 

For every fc > 1, the k*^ construction step is Wjt+3<i(t'fc+2. v^.), which appends one vertex 
and two edges to the graph, and forms a nested quadrilateral Qk = Vk+3Vk+2Vk+iVk- We 
denote the four edges of Qk as {vk+3,Vk+2) = Sfc.i, (ufe+2, ffe-i-i) = Sk,2, ivk+i,Vk) = Sk.3 
and (ufc, Ufc+a) = Sfe,4, and two diagonals (wfe+3, Wfe+i) = fk+i, {vk+2,Vk) = fk- Notice that 
Qk shares two edges with Qk-i - Sfe,2 = Sfe-1,1, Sfc,3 — Sfc-i,2. The forward realization type 
a is assigned such that Vk+3 and Vk+i lies on different side of fk- 

Clearly, G is 1-path with low Cayley complexity. So we use the QIM algorithm intro- 
duced in Section 4.1.1 to compute <1?/(G, T). We start from the last extreme edge, and 
repetitively map l{fk) to get intervals for l{fk-\), until we obtain the intervals for 



V2 




Figure C.l. For Observation C.l: a 1-dof tree-decomposable graph with exponential Cayley size 
for base non-edge /i = {vi^vs). The graph is a series of nested quadrilaterals and each cluster 

is an edge. 

We choose Tsuch that for Qi, r(si,i) = l{v4,v^) = 8, l{si^2) = 1(^3, V2) = 8.1, ^(si.s) = 
l{v2,vi) = 7.9, r(si,4) = l{v4,vi) = 1. The ellipse C relating the two diagonals of Qi is 
shown in Figure C.2. Only the upper half of the curve (shown in solid line) corresponds 
to realizations with forward realization type a. 

For each fc > 1, we assign l{sk,i) and l{sk,4) by the following: (1) Observe the ellipse 
C of Qk-i, as shown in Figure C.2. Denote the leftmost point on C by Pmin(/fc-i); the 
rightmost point by Pmax(/fe-i), the topmost point by Pmax(/fe)- Let I2 be the length of fk 
at Pmax(/fc), and h be the larger one of the lengths of fk at Pmin(/fe-i) and Pmax(/fe-i)- 
The interval [hjh] is attainable by l{fk) in Gf{k — 1). (2) Assign l{skj) and l{sk.4) such 
that l{8k,i) - l{skA) = r™(/fc) = '(1 + e)h, l{sk,i) + l{sk,4) = ^'"^''(/fe) = (1 - e)^2, where 
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e is a positive value small enough so that l{sk,i) and l{sk,4) have positive solutions. In 
this way, l{fk) is restricted to an interval [l'^^'^{fk)J'^'^^ifk)] slightly tighter than [lijh]- 
For example, we want to assign I for the 2"^^ realization step ^5 < {v4,V3). As shown 
in Figure C.2, in Qi, h = 7.9 + 1 = 8.9. For p^inifk-i), Kfi) = 8 - 1 = 7, the 
corresponding l{f2) ~ 8.36. For Pmax(/fe-i), Kfi) = 8 + 1 = 9, the corresponding /(/2) ~ 
7.42. Therefore h = nmx[8.36, 7._42] = 8.36. So let Z™(/2) = /(w5,W4) - Z(w5,W3) = 
(1 + 10-5);^ « 8.364, r''''lf2) = 1{V5,V4) + Iiv5,v3) = (1 - 10-^)12 ~ 8.900. We assign 
l{v5,Vi) « 8.632, 1{V5,V3) « 0.268. The two new extreme linkages ((7/(2),?™) and 
((5/(2),/™^^) each has two realizations, and each of these realizations creates a new 
endpoint in ^f{G,l): (G/(2), P^'^'^) corresponds to realizations in Figure C.2 (b) and 
(d) , ((7/(2), r™'") corresponds to reahzations in Figure C.2 (a) and (c) (point b, d, a 
and c in the left graph respectively). Since lb{fi) ~ 7.00, ld{fi) ~ 8.52, Za(/i) ~ 7.49, 
Icifi) ~ 7.51, $/(G,0 contains two intervals h = [7.00,7.49] and h = [7.51,8.52], 
corresponding to two different reverse realization types. 



/2 = {V4,V2) 





a- 
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"fl = {vi,V3U^ 





V4 (d) 



Figure C.2. Ellipse C for quadrilateral Qi = «4«3«2«i. Length of extreme edge {v2,V4,) is re- 
stricted by the 2""* realization step, and of l{vi,V3) has 2 intervals. 

Table C.l shows I for the subsequent construction steps, computed by the procedure 
described above. Figure C.3 shows $/(G/(6),r). The single interval of l{v5,V3) maps 
to 2 intervals for ?(u2,W4): h = [8.36,8.48], h = [8.49,8.74], and 4 intervals for the 
base non-edge Hvi.vs): In = [7.000,7.008], hi = [7.010,7.121], = [8.039,8.391], 
/12 = [8.403,8.524]. 

In general, the fc*'' realization step (fc > 1) produces one interval for l{fk) which 
maps to 2 intervals for l{fk-i), 4 intervals for l{fk-2), and finally 2'"'"^ intervals for 
fl = (I'liVs)- Notice that there is no overlapping since the curve is monotonic in each 
interval. 
□ 
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number of intervals for l{vi, 113) after step k 


2 


8.632 


8 


8.1 


0.268 


2 


3 


8.306 


8.632 


8 


0.062 


4 


4 


8.044 


8.306 


8.632 


0.017 


8 


5 


8.645 


8.044 


8.306 


0.004 


16 


6 


8.310 


8.645 


8.044 


0.001 


32 


7 


8.045 


8.310 


8.645 


0.0003 


64 


8 


8.645 


8.045 


8.310 


0.00006 


128 


9 


8.310 


8.645 


8.045 


0.00001 


256 


10 


8.045 


8.310 


8.645 


0.000004 


512 



Table C.l. Edge lengths of quadrilateral Qk for construction step 2 to 10. 
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Figure C.3. Ellipse C for quadrilaterals (a) Q2 ~ V5V4V3V2 and (b) Qi = after 
realized. The Cayley configuration space over {vi, V3) is divided into 4 intervals In, I21, I22, I 
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